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Abstract 

In order to find analytically the travelling waves of partially integrable au- 
tonomous nonlinear partial differential equations, many methods have been pro- 
posed over the ages: "projective Riccati method", "tanh-method" , "exponential 
method", "Jacobi expansion method", "new ...", etc. The common default to all 
these "truncation methods" is to only provide some solutions, not all of them. By 
implementing three classical results of Briot, Bouquet and Poincare, we present an 
algorithm able to provide in closed form all those travelling waves which are elliptic 
or degenerate elliptic, i.e. rational in one exponential or rational. Our examples 
here include the Kuramoto-Sivashinsky equation and the cubic and quintic complex 
Ginzburg-Landau equations. 
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1 Introduction 

In physical language, given an autonomous nonlinear partial differential equation (PDE) 
E(v , v x , Vt, ■ ■ ■ , v nXjmt ) = 0, one calls travelling wave any solution of the reduction 
(v,x,t) — > £ = x — ct,v(x,t) = w(£), with c constant, and solitary wave a travel- 
ling wave which obeys the decay condition that u(£) should have constant limits, possibly 
different, as £ — > ±oo. The problem addressed here is to find all the travelling waves of 
a given PDE, in closed form, or more generally all the solutions of similarity reductions 
of a given PDE to some autonomous ODE. Two cases naturally arise. If the PDE is 
"integrable", for instance in the sense of the inverse spectral transform [Tj, there exist 
powerful methods to solve this problem completely; we will discard this case. If the PDE 
is partially integrable (i.e. fails at least one of the integrability criteria), there exist several 
nonperturbative methods (mostly "truncations" [28] ) able to find some travelling waves 
but generically unable to find all of them. 

In a more mathematical language, let us denote N the differential order of the au- 
tonomous ordinary differential equation (ODE) 

E(u,u',...,u^)=0, ' = A (1) 

which is the travelling wave reduction of the considered PDE. If the ODE has the Painleve 
property [TT] (i.e. no "bad" singularities), the problem is to find its general solution 
(depending on N arbitrary constants) in closed form; we discard this case. If the ODE fails 
the Painleve test ("partially integrable" case), by definition it admits some movable critical 
singularities (movable: which depends on the initial conditions; critical: the solution is 
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multi- valued around it). Let us define its general analytic solution as the M-parameter 
particular solution without movable critical singularities, with M maximal and of course 
M < N. The problem is then to find this general analytic solution in closed form. This 
is equivalent to find the largest order subequation of the given ODE with the Painleve 
property (subequation: the ODE is a differential consequence of it). 

The present article, which is self-contained, gives a synthetic presentation of a method 
[2"B1 [TT] able to find this general analytic solution (i.e. all the single-valued travelling 
waves) when it is elliptic or degenerate of elliptic. This includes the physically relevant 
solitary waves. 

There are many equivalent definitions for elliptic functions. From the point of view 
of nonlinear differential equations adopted here, the suitable definition for an elliptic 
function /(£) is any rational function of p'{£), 

m) = Rm),m), (2) 

in which p is the function of Weierstrass, defined by the first order ODE 

p' 2 = V - g 2 p -g 3 = 4(p(0 - ei )(p(0 - e 2 )(p(0 - e 3 ), p" = 6p 2 - |. (3) 

Elliptic functions are doubly periodic and have two successive degeneracies, 

• when one root ej is double (<?f — 27 g 2 = 0), degeneracy to simply periodic functions 
(i.e. rational functions of one exponential e kx ) according to 



Vx, d : p(x, 3d 2 , -d 3 ) = 2d - 




• when the root tj is triple (g 2 — #3 = 0), degeneracy to rational functions of £. 

In section [21 we handle a textbook example of an integrable equation in order to make 
several points. In section [3j we make precise the notion of general analytic solution of a 
partially integrable ODE. In section HJ we shortly review the existing methods and point 
out their weaknesses. In section [5l we present the algorithmic method, then we process a 
few selected examples. Finally, section [7] examines possible extensions of the method. 



2 Tutorial integrable example 

Consider the Korteweg-de Vries (KdV) equation 

6 

bvt + v xxx vv x = 0, (a, b) constant, (5) 

and its travelling wave reduction 

oh ... 6 . 

v = c + u, t = x — ct, u uu = 0. (6) 

6 a 

2.1 Mathematical part 

If one a priori knows nothing on ([6]), some skill is needed to find the successive first 
integrals g 2 (arising from the conservative form) 

3 

u" u 2 + ag 2 = 0, (7) 
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and gs (integrating factor u') 



^u' 2 - —v? + ag 2 u + 2a 2 g 3 = 0. (8) 



The general solution is then elliptic (doubly periodic) 

u(0 =2ap(£- £ 0} g 2 ,g 3 ), (9) 
in which the three constants of integration ^o,fl'2,5 , 3 are allowed to be complex. 

2.2 Physical part 

One then prescribes some behaviour at £ — > ±oo, e.g. 

lim u = B_, lim u = B + , B± real, (10) 

which for B_ = B + = B makes the doubly periodic function degenerate to simply periodic 
(trigonometric) 

g 2 = 3{{B - K)/a)\ g 3 = -{{B - k)/o)\ k = -abc/6. 
u = -^c + 2a (h> - coth 2 k{i - &)) , k 2 - 3{B ~ K) (U) 



6 V3 v V 2a 

in which the constants of integration are k 2 (real) and £o (complex). 

Finally one requires this solution u to be bounded on the real axis, by moving all 
singularities outside 1Z 

£o = £i+"r/(2A;), 6 e K, cothfc(f -£o) = tanhA:(£ - &), (12) 
to obtain the physical, bell-shaped, solitary wave on the background B, 

u = B- 2ak 2 sech 2 (kx - AkH/b + 6k(B/a)t/b - jfefr) . (13) 

2.3 Lessons from this tutorial example 

Although it is completely straightforward, this example nevertheless teaches us several 
lessons. 

• Skill was needed to find the first integrals, we want to avoid that. 

• If possible, it is better to first find the mathematical solution, because all the physical 
solutions follow. In the case of the nonlinear Schrodinger equation for instance, 
which admits two (one bright, one dark) solitary waves depending on the focusing 
or defocusing situation, the bright wave is found as above, and the dark one arises 
from a factorisation of the first order subequation analogous to ([8]). 

• Elliptic functions seem privileged. We will come back to that later. 
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3 Partially integrable case 

By definition, partially integrable PDEs fail the Painleve test. Our specific examples will 
be the following reductions: 

• the similarity reduction of the one-dimensional cubic complex Ginzburg-Landau 
equation (CGL3) 

%A t + pA xx + q\A\ 2 A-i 1 A = 0, pg 7 ^ 0, (A, p,q)eC, 7 G K (14) 
A(x, t) = ^M{i)e^- Ut + ^(0) , £ = x -ct, (15) 
(a generic equation for slowly varying amplitudes, see the review [25J); 

• the travelling wave reduction of the Kuramoto-Sivashinsky equation (KS) 

v t + w xxxx + bv xxx + fxv xx + vv x = 0, v ^ 0, (16) 

u 2 

v(x, t) =c + tt(f ), £ = x-ct, vv!" + bu" + fm' + — + A = 0, v ^ 0, (17) 

in which A is an integration constant (flame on a vertical wall, phase equation of 
CGL3, see (22] for some physical background); 

• the similarity reduction of the one-dimensional quintic complex Ginzburg-Landau 
equation (CGL5) 



%A t + pA xx + q\A\ 2 A + r\ A\ 4 A - i^A = 0, pr ^ 0, 9f(p/r) ^ 0, (18) 
(A,p,q,r) eC, ieTZ. 

A(x,t) = y/M(0e i{ ~ ut + ^\ £ = x-ct. (19) 

3.1 General analytic solution of a partially integrable ODE 

Given an iV-th order ODE (DO) which is partially integrable and therefore displays movable 
multi-valuedness, one first needs to count the number of integration constants which 
correspond to this multi-valuedness. Let us take an example. 

3.2 Local representation of the general analytic solution 

The KS PDE ( JT6l) admits the travelling wave reduction (fTTl) . It has a chaotic behavior [22], 
and it depends on two dimensionless parameters, b 2 /(fiu) and uA/fi 3 . The assumption of 
a singular algebraic behaviour 

£ ^ £0 : it ~ u oX P , X = Z ~ Co, u ^ 0, p i M, (20) 
generates for (p, it Q ) the system 

p - 3 = 2p, up(p - l)(p - 2)«o + (l/2)«g = 0, u ^ 0, (21) 

which admits the unique solution p = —3, uo = 120z/, with the common value of the two 
powers: q = p — 3 = 2p = —6. Near this triple pole, the linearized equation has the 
indicial equation 

\\m X - j - q {vdl + u QX p )x j+P (22) 
= v(j - 3)(j - 4)0' - 5) + 1201/ = i/O' + 1)0 2 - 13j + 60) (23) 

= ^ + 1) 6 - H±i^ _ Hzi^ = o, (24) 



i.e. two of the three Fuchs indices are irrational complex numbers. Computing the next 
terms after the dominant one yields the Laurent series 

(o) im 3 ir , 2 15(16/xi/ - b 2 ) _! 13(4/iz/ - 6 2 )6 ^. . 

which only depends on one integration constant. Laurent series such as (1231) are conver- 
gent, as proven by Chazy [3]. By perturbing this solution [3, pages 338, 384] [7J, one 
obtains the three-parameter solution 

«(£-£o,£C+,£C_) = M (0) 

+^X" 3 [c-iX~ 1 Regular(x) 

+ c + x (13+iVYT)/2 Regular( X ) 

+ c_x (13 - iv/7T)/2 Regular(x)] + C(£ 2 )}, (26) 

which is a local representation valid for £ — > £o and £ — > 0, depending on the three 
arbitrary constants £ Q , £c + ,£C_ (indeed, as shown by Poincar43 and Darboux [13], the 
coefficient of ec_i in ([26]) is the derivative of with respect to its arbitrary constant £o ; 
therefore c_i only represents a perturbation of £o an d it can be set to 0). The only way 
for (1261) to become single-valued is to cancel the contribution of the two irrational Fuchs 
indices, thus restricting (|23|) to (1231) . 

The general analytic solution as defined in section [TJ depends on 3 — 2 = 1 integration 
constant, it is locally represented by the Laurent series (1231) and the problem addressed 
is to obtain a closed form expression for u^ '. 



4 Existing methods (sufficient) 

Given a nonlinear ODE such as (pQ), all the methods able to yield single- valued solutions in 
closed form necessarily make use of the structure of movable singularities. Many "new" , 
"extended", etc, methods claim not to do so and to be original, but this is not true: no 
method can avoid the analytic structure of singularities. 

Moreover, in order to be nonperturbative, i.e. to only require a finite number of steps, 
the integration methods have better use an elementary "unit of integration" which has 
singularities (such as tanh or sech) rather than an entire function (such as exp). 

When nothing is known about the ODE (TjQ), most methods follow the pioneering 
work of Weiss, Tabor and Carnevale [28J. The restriction of their method to ODEs 
was certainly known to the classical authors involved with explicit integration (Painleve, 
Gambier, Chazy, Bureau), although we could not find an explicit, convincing reference. 

All these methods proceed as follows: 

1. Assume a given class of expressions for the general analytic solution, 

2. Check whether there are indeed solutions in that class. 

Various names are used to describe these methods: truncation method, tanh method, 
Jacobi expansion method, .... 

Typical classes of assumed expressions are: 

• polynomials in tanh (Weiss et al. [28] and followers [5]), 

• polynomials in p, p' pH EH] , 

• polynomials in tanh and sech (TTJ [HI [24"] . 

We will call all these methods sufficient because, by construction, they cannot find 
any solution outside the given class of expressions. 

1 With his "equation aux variations" , see Chazy [31 page 338 footnote 1] . 
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4.1 Examples of sufficient methods 

The KS ODE (fTTj) has one family of movable triple poles, and the Weierstrass function p 
one family of double poles. Therefore it is consistent to assume 



u = c p + Cip + c 2 , c ^ 0, 



u 



E(u) = vu" + bu" + fiu + — + A = 0, v ^ 0, 
P =6p -y, p =4p -g 2 p-93- 
The elimination of p" and p' 2 transforms the lhs of ( |27l) into 

j=0 fc=0 

and one has to solve the six determining equations 

f £ 3 ,o = c (120z/ + 2co) =0, 
Ei i = 126c + c cx + 12z/ci = 0, 



(27) 



2.0 



6/iCo + 66ci 



0. 



#o,i = c c 2 + /iCi = 0, 
#i,o = cic 2 - -5- 2 Co - I8ug 2 c = 0, 

^o,o =A - 12ug 3 c + - (c 2 - 6^ 2 ci - /i# 2 c - flfccg) = 0. 

As always in these truncations [28J, these equations mttsi be solved by decreasing value 
of the singularity degree 3j + 4k, for the successive unknowns c , ci, c 2 , <? 2 , #3- The result 

is iiiisa 



it = — 60z/p' — 156p 



6/i 

4^' 



#2 



/i- 



12z/ 2 



#3 



1080z/ 3 ' 



(2? 



16/iz/, 



and this would be the general analytic solution if b 2 — 16/iz/ were unconstrained. Removing 
this constraint is still an open problem. 

Instead of considering the Weierstrass equation, one may consider the Riccati equation 
with constant coefficients 



r' + r 2 + !=0, S 



2 



constant G C, 



(29) 



whose singularities are one family of movable simple poles, and make the similar consistent 
assumption [3J [8] 



i=o 



CjT 



-J-P 



Co ^ 0, p 



(30) 



with Cj constants to be determined. The elimination of r' from f[2"9"j) transforms the lhs of 
(|2T|) into a similar polynomial 



3-1 



0, g 



-6. 



(31) 



j=0 



This set of determining equations Ej = admits six solutions, listed in Table [TJ all 
represented by 



u = nOvr* - 156r 2 + I ^-u - 30uk z - - -— ) r 

I 19 4 x 19v] 



60 „ 15b 2 



5 Tl9 13£» 3 7 rib k , k tJ 

+ -bk 2 + " , r = - tanh -(£ - £ ). 32 

2 32 x 19z/ 2 4xl9z/' 2 2^ 



The solitary waves 6 = were found by Kuramoto and Tsuzuki [2T], and the three other 
values of b 2 /(/iu) were added by Kudryashov [19] . 
Currently, no other solution is known to fl27l) . 



Table 1: The six trigonometric solutions of KS, Eq. (ITT)) , with the notation k 2 = —2S. 
They all have the form ()32l) . The last line is a degeneracy of the elliptic solution ( l28l) . 



67(f) 


iM//i 3 


isk 2 /fi 





-4950/19 3 , 450/19 3 


11/19, -1/19 


144/47 


-1800/47 3 


1/47 


256/73 


-4050/73 3 


1/73 


16 


-18, -8 


1, -1 



It is easy to build an example making truncation methods fail. Knowing that rational 
functions are usually not assumed classes, such an example is for instance 



u = - — TT77 ~ = ou ^ s i c i e usua l classes, (33) 



tanh(£ - go) 

2 + tanh 2 (e-eo) 
2u' 2 + (24m 2 - 3)u + 72m 4 - 17m 2 + 1 = 0. (34) 



5 Subequation method (necessary) 
5.1 The class of elliptic functions 

All solutions found in previous examples are elliptic or degenerate. This quite frequent 
result for solutions of autonomous nonlinear ODEs has a simple explanation, dating back 
to Lazarus Fuchs. Indeed, this is a classical result (L. Fuchs, Poincare, Painleve) that 
the only autonomous first order algebraic ODEs having a single-valued general solution 
are algebraic transforms of either the Weierstrass elliptic equation, see ([2]), or the Riccati 
equation ( 129)) . These first order ODEs are therefore the elementary building blocks with 
which to build solutions to autonomous ODEs of higher order. 

Conversely, given an iV-th order autonomous algebraic ODE, is it possible to find 
all solutions in this class of elliptic or degenerate elliptic functions? The present section 
brings a positive answer. 



5.2 Two results of Briot and Bouquet 

Given an elliptic function, its elliptic order is defined as the number of poles in a period 
parallelogram, counting multiplicity of course. It is equal to the number of zeros. 
We will need two classical theorems. 
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Theorem theorem XVII p. 277]. Given two elliptic functions u, v with the same 
periods of respective elliptic orders m, n, they are linked by an algebraic equation 

m n 

F{u,v) = Y,T l HkU i v k = 0, (35) 

fc=0 j=0 

with deg(F, u) = order (v), deg(F,v) = order(w). If in particular v is the derivative of u, 
the first order ODE obeyed by u takes the precise form 

m 2m— 2k 

F(u,u') = J2 J2 a j,kU 3 u k = 0, a , m ^ 0. (36) 

k=0 j=0 

Theorem (Briot and Bouquet, Poincare, Painleve). If a first order m-th degree au- 
tonomous ODE has the Painleve property, 

• it must have the form (|36|) . 

• its general solution is either elliptic (two periods) or rational in one exponential e kx 
(one period) or rational in x (no period) (successive degeneracies g\ — 27 g\ = 0, 
then g 2 = in p' 2 = 4p 3 - g 2 p - g 3 ). 

Remark. Equation (|36j) is invariant under an arbitrary homographic transformation 
having constant coefficients, this is another useful feature of elliptic equations. 



5.3 General method to find all elliptic solutions 

Consider an iV-th order autonomous algebraic ODE (pQ) admitting at least one Laurent 
series 

+oo 

u = X p X = £ - fo- (37) 

3=0 

If its general analytic solution, as defined section [TJ is elliptic or degenerate elliptic, 
there exists an algorithm to find it in closed form, which we now presen t. It only requires 
four ingredients: the two theorems of Briot and Bouquet stated section 15.21 the existence 
of at least one Laurent series, an algorithm of Poincare to be presented soon. Its input 
and output follows. 

Input: an iV-th order (N > 2) any degree autonomous algebraic ODE admitting a 
Laurent series. 

Output: all its elliptic or degenerate elliptic solutions in closed form. 
The successive steps are (231 IH] : 

1. Find the analytic structure of singularities (e.g., 4 families of simple poles, 2 of 
double poles). Deduce the elliptic orders m, n of u, v! . 

2. Compute slightly more than [m + l) 2 terms in the Laurent series. 

3. Define the first order m-th degree subequation F(u, u') = (it contains at most 
(m + l) 2 coefficients 

m 2m— 2k 

F(u,u') = Y^ Yl a j,kU 3 u' k = 0, a 0jTO ^ 0. (38) 

k=0 j=0 
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4. Require each Laurent series fl37j) to obey F(u,u') = 0, 

F = X m{p - 1] F rf + 0(x J+1 ?j , Vj 

and solve this linear overdetermined system for aj^ 

5. Integrate each resulting ODE F(u,u') = 0. 
Several remarks are in order. 

1. The fourth step generates a linear, infinitely overdetermined, system of equations 
Fj = for the unknown coefficients cij t k- It is quite an easy task to solve such a 
system, and this is the key advantage of the present algorithm. 

2. In the fifth step, two cases arise, depending on the genus of the algebraic curve 
F — 0. If the genus is one (nondegenerate elliptic), there exists an algorithm of 
Poincare which builds explicitly the rational function R, Eq. (12"1). This algorithm 
has been implemented by Mark van Hoeij [15] as the package "algcurves" of the 
computer algebra language Maple. If the genus is zero (degenerate elliptic), finding 
u as a rational function of one exponential e fc ^ - ^ ^ or as a rational function of £ — £o 
is classical, and also implemented in "algcurves". 

3. In the third and fourth steps, the requirement can be weakened, at the price of 
finding less that the general analytic solution. Consider for instance the artificial 
ODE (which actually occurs in the travelling wave of the modified KdV equation), 

E{u) = some differential consequence of a 2 u' 2 — (u 2 + b) 2 + c = 0, (40) 
which admits two Laurent series 

u = ±ax~ 1 + .... (41) 

With the assumption m = 2 in step 3 and "Require each Laurent series" in step 4, 
one will find the Jacobi solution. But, with the weaker assumption m — 1 in step 
3 and "Require one [not two] Laurent series" in step 4, one will find the constraint 
c = in E(u). The interested reader can practice on the ODE admitting, for 
instance, the rational solution 

u = ci(x — xi)^ 1 + c 2 (x — X2)~ 2 + (x — xs)~ 3 , Cj fixed, Xj movable, (42) 
for which the assumptions can be m = 6, 5, 4, 3. 



F, = 0. (39) 



6 Examples 

6.1 Tutorial integrable example: KdV 

The ODE ([6]) presents one movable double pole u ~ 2ax~ 2 , hence the respective elliptic 
orders 2 and 3 for u and u' . Nine terms (in fact five because of parity) will prove sufficient 
in the computation of the unique Laurent series, 

u = 2a X - 2 + U AX 2 + U eX 4 + ^X 6 + • • • , (43) 
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in which U4 and Uq are arbitrary constants. Assuming in step 3 



F = u' + a Q iu' + ai^uu' + a ,o + «i,o« + «2,oW 2 + a3,oU 3 , a , 2 = 1, (44) 
one generates in step 4 the linear overdetermined system ( 1391 ). 

F = 16a 2 a , 2 + 8a 3 a3,o = 0, 
iq = -8a 2 a hl = 0, 
F 2 = 4a 2 a 2>0 = 0, 
F 3 = -4aa ,i = 0, 

F 4 = 2aai. — 16aao.2^4 + 12a 2 a 3)0 f/4 = 0, 

F 6 = a 0i0 + 4aa 2 ,o^4 - 32aa 0i2 f/ 6 + 12a 2 a 3i0 L/" 6 = 0, 



whose unique solution is 

u' 2 - (2/a)u 3 + 2OC/4W + 56aU 6 = 0. (46) 

The two arbitrary constants C/4, correspond to the two first integrals of ([6]), but no skill 
is required to find them, the process is systematic. 
In step 5, the Maple commands would look like 

with(algcurves) ; load the package 

genus (eq46 ,u,uprime) ; will answer "1" 

Weierstrassf orm (eq46 , u , uprime , wp , wpprime) ; 

the last command yielding the four formulae of the birational transformation between 
(u,u') and (p, p'), one of them being precisely Q). 



6.2 Partially integrable example, KS 

The Laurent series of ( |T71) is ( 1251) . 

In the first step, the unique family of movable triple poles yields elliptic orders 3 and 4 
for respectively u and u'. With the normalization 00,3 = 1, the subequation F(u, u') = 
contains ten coefficients. In step 4, these ten coefficients are first determined by solving 
the Cramer system of ten equations Fj = 0,j = : 6,8,9,12. The remaining infinitely 
overdetermined nonlinear system for (u, b, fi, A) contains as greatest common divisor (gcd) 
b 2 — 16fiu which defines a first solution 

— = 16, u s = u H -, 

jju 32u 2 

t +—) 2 = °- (47) 
1024i/ 4 3 J K J 

After division by this factor, the remaining nonlinear system for (u, b, /i, A) with b 2 — 
Wfiu 7^ admits exactly four solutions (stopping the se ries at j = 16 is enough to obtain 
the result), identical to those listed in Table [Tl Section [4.1} each solution defining a first 
order third degree subequation, 

6 = 0, 

/ 18(V 2 \ 2 /, 360y\ 9 / 2 30>, 30VV n , Asl , 
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, u s = u- 




, n /3 9 / , 30u . 30V i 

b = 0, u' 3 + — ( M 2 + -£u' + ] = U. ( 19) 

6 2 144 



^ 47 > ^ = u - 7773 > I M ' + — M * ) + —. u s = 0, (50) 
b 2 _ 256 

~ - "73"' -» - 2048i/ 2 




/ & \ 9 f 2 5&3 5&2 A n 

w H m s H u. H ;ru s H w = 0. (51) 

2z/ 7 40i/\ s 1024^ 2 s 128f J V ; 

In order to integrate the two sets of sub equations (147] ) , P5|) -( |5T|) . one must first 
compute their genus, which is one for (J47J, and zero for (148]) - (|51|K Therefore (1471) has the 
elliptic general solution ( 128]) . 

As to the general solution of the four others ( 148]) -( l5Tl) . this is the third degree poly- 
nomial ( |32l) in tanh&;(£ — £o)/2 already found by the one-family truncation method. 

Remark. Canceling the gcd b 2 — 16/ii/ is equivalent to cancel the residue in the Laurent 
series (125]) . therefore KS admits no elliptic solution other than (1281 [16J. Indeed, this is a 
property of elliptic functions that, inside a period parallelogram, the sum of the residues is 
necessarily zero. Conversely, by implementing the condition that the sum of the residues of 
any differential polynomial of u vanishes, Hone [16J has found an elegant method to isolate 
all the nondegenerate elliptic (genus one) solutions, proving in particular that CGL3 has 
no travelling wave which is nondegenerate elliptic. This method cannot however also yield 
the degenerate elliptic solutions. 

6.3 Partially integrable example, CGL5 



As example of a new solution found by the method of section 15.3] one can quote the 
unique elliptic (nondegenerate, i.e. genus one) travelling wave ( TT9|) of the CGL5 equation 

M V(Q - 92 , cs,( 3g 2 \ 1/2 

notation s r — isi = p , g r = s r 7 + SjU;, = $S(r/p), ^2) 
constraints q — 0, %t{r/p) = 0, csi = 0, s r u — + (cs r /2) 2 = 0, 



values of c ,g 2 ,g3 ■ c = — , g 2 = -— , #3 = 0. 



„2 _ 

" 3e,-' aA 27'' 



This solution, obtained by combining the present method with the conditions on residues 
|16j . exists at the price of five real constraints among the coefficients p, q, r, 7, u, c of the 
differential system for (M, if') . 

In this example, care should be taken that the two elliptic subequations, 



3 4 e;M' 4 - M 2 (3e;M 2 - Ag r f = 0, 

9^ 2 -12^-^ 2 = 0, ^ = ^'-f : , (53) 

respectively of a canonical Briot-Bouquet type and of a Jacobi type, involve elliptic func- 
tions with different first arguments, see p(2£) term in 



6.4 Results for various partially integrable PDEs 

For the various PDEs considered in the text, and recalled below for convenience, 
KdV : u'" - (6/a)uu' = 0, 
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KS : vu"' + bu" + /m' + u 2 /2 + A = 0, v ^ 0, 
CGL3 : iA* + pA xx + q\A\ 2 A - i 1 A = 0, pq ^ 0, 3(g/p) ^ 0, 
CGL5 : iA t +pA xx + q\A\ 2 A + r\A\ A A-i~?A = 0, pr ^ 0, 9f(p/r) ^ 0, 
2CGL3 : S^At = rA ±T ^v4 ± + (l + ^a)^A ± -(l + z/3)(|v4 ± | 2 + 7 |A T | 2 )A ± , 
r, v, a, f3, 7 real parameters, 

Table [2] collects the results produced by the method for the various PDEs in the text (we 
have added in last line an interesting system, not yet processed). In this table, nPq is 
short for "n families of q-th order poles" , (m + 1) 2 — 1 is the maximal number of coefficients 
in the subequation (I36I) . and the column shows the true number of such coefficients. 

The column displays the minimal number of terms to be computed in the Laurent 
series ( 1371) . 

Table 2: Summary of results. The solutions are labelled "ell, trig, rat" for, respectively: 
elliptic (genus one), rational in one exponential, rational. 



PDE 


poles 


m 


(m+ l) 2 - 1 


a j,k 


Uk 


solutions 


Ref 


KdV 


1P2 


2 


8 


6 


8 


1 ell 




KS 


1P3 


3 


15 


10 


16 


1 ell + 4 trig 


m\ 


CGL3 


2P2 


4 


24 


18 


24 


6 trig + 1 rat 


[23,116] 


CGL5 


4P1 


4 


24 


24 




1 ell + n trig 


[27j 


2 CGL3 


2P2 


4 


24 


18 




To be done 


ma 



7 Possible extensions of the method 

It should first be noted that the output of the method is all solutions which are elliptic 
functions of £ — £o or rational in one exponential e k ^~^ or rational in £ — £ 0) i n which £ 
is the arbitrary constant arising from the first order subequation. 
Having this in mind, one can think of several possible extensions. 



7.1 Extension to autonomous discrete equations 

If the independent variable £ is discrete (£ = nh, with n integer and h a given stepsize), 
can the method be extended so as to yield all elliptic or degenerate elliptic solutions? The 
question is still open. 

For instance, the discrete nonlinear equation of the Schrodinger type 



u(x + h, t) + u(x — h, t) — 2u(x) \u\ 2 u 

iu t + p— — y — — — + q — 7 = 0, pqv ^ 0, 

h 2 q i + u(qh 2 /p)\u\ 2 



(54) 



admits various elliptic solitary waves [T8j [6] but no proof exists that these are the only 
ones. 

Note that it is sufficient to look for rational functions of £>(£),£>'(£) since, as shown 
by Briot and Bouquet [2], any elliptic function of argument £ + a can be expressed as a 
rational function of p'(£)- 



13 



7.2 Extension to nonautonomous differential equations 



Consider for instance the problem (already solved, but this is just an example) of finding 
by the above method all the rational solutions of the second Painleve equation P2 

u" = 2u 3 + xu + a, (55) 

i.e. the sequence, for a integer, 

1 . . 

u = =f— , a = ±1, (56) 

x 

2(x 3 — 2) 

" = *^ry a = ±2 ' (57) 

3x 2 (x 6 + 8x 3 + 160) lo /eo . 

u = =p— — — — -, a = ±3, .... (58) 

(x 3 + 4)(x 6 + 20x 3 - 80) V ; 

The integration constant £o of previous sections is no more arbitrary. 

In order to retrieve these rational solutions by the method of section 15.31 one must 

• ignore the structure of movable singularities of ( 1551) ; 

• consider instead the simplest differential consequence of fl55l) which is autonomous, 
i.e. (elimination of x), 

uu'" - u'u" - Au 3 u' - u 2 + au' = 0; (59) 

• guided by (156|) - (|58|) . manage to find the (probably a 2 in number) families of movable 
simple poles in (I59|) (it is for the moment not completely clear how this result should 
come out); 

• apply the method to fl59|) with m = a 2 and find a unique solution. 

We have not yet done it but these are at least the guidelines to be followed. 
Remark. The symmetries of P2 allow to considerably simplify the above by the change 
(u,x) — > (U,X),u = x 2 U,X = x 3 , thus lowering m to \a\. 



7.3 Extension to Painleve solutions 

After the elliptic function, the next elementary functions defined by differential equations 
are the six Painleve functions, which all obey a nonautonomous nonlinear second order 
ODE. Hen ce the question: given an iV-th order ODE (OQ), iV > 3, can the method of 
section IBTBI be extended so as to yield all the solutions which are algebraic transforms of a 
given Painleve equation? Let us take an example to illustrate the difficulties to overcome. 

Consider the third order autonomous ODE built by elimination of the constant K\ in 
(this example is taken from the Lorenz model) 

d 2 x dx x 3 /8 Ki 9 ,\ . , 

+ 2— + — + - - -i e - 2 ' x = 0. (60) 



dt 2 dt 2 V9 2 
This ODE happens to be integrable with a particular second Painleve function with a = 

x = ae- 2t / 3 X } T = V^e" 2 '/ 3 , K x = ^a 3 , (61) 
2 8 

= 2X 3 + TX, (62) 
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The problem is to uncover the subequation fl62|) for X(T) knowing only the third order 
autonomous ODE for x(t) and its Laurent series. The necessity to take account of a 
possible change of the independent variable t — > T introduces a countabl^ number of 
arbitrary coefficients (those of the Taylor series oit^T near the movable singularity of 
x(t)) and thus makes it difficult to use the information provided by the Laurent series. 

In fact , th e main privilege of the elliptic functions, already mentioned as a remark 
in section 15.21 is the autonomous nature of their ODE and the invariance of this ODE 
under homographic transformations, all features which do not exist any more for Painleve 
functions. 



8 Conclusion and perspectives 

The main features of the method presented here are as follows: 

• It provides all the elliptic and degenerate elliptic (i.e. rational in a single exponential 
or rational) solutions of a given autonomous nonlinear algebraic ODE; 

• It includes all "truncation" , "extended" , "new" , etc methods; 

• It makes all these methods obsolete; 

• Its main part (solving system ( 1391) ) is a linear algebraic problem. 

Of course, integrating the resulting first order ODE may consume much more time. 

There exists another theory in which meromorphic functions play a central role, this 
is Nevanlinna theory. For a short comparison between Nevanlinna theory and the present 
method, the interested reader can refer to [12] . 

Current and future work includes: 

• Revisit physically interesting PDEs with i nsuffi cient solutions, such as the system 
of two coupled CGL3 mentioned in section I6.4L 

• Find general analytic solutions which are not elliptic, for instance Lame functions. 

Acknowledgements 

Partial financial support has been provided by the Research Grants Council contracts 
HKU 7123/05E, HKU 7184/04E and HKU 7038/07P. RC acknowledges the partial sup- 
port of the PROCORE - France/Hong Kong joint research grant F-HK29/05T. 

References 

[1] M. J. Ablowitz and P. A. Clarkson, Solitons, nonlinear evolution equations and in- 
verse scattering (Cambridge University Press, Cambridge, 1991). 

[2] C. Briot et J.-C. Bouquet, Theorie des fonctions elliptiques (Mallet-Bachelier, Paris, 
1859). http : //gallica. bnf . f r/document?0=N099571 

[3] J. Chazy, Sur les equations differentielles du troisieme ordre et d'ordre superieur dont 
l'integrale generale a ses points critiques fixes, These, Paris (1910); Acta Math. 34 
(1911) 317-385. 

2 This number can be made finite by requiring the change (x, t) — > (X, T) to introduce only exponential 
functions as in (|6f p . so as to make the transformed ODE polynomial in the new independent variable T, 
before application of the method. 



15 



[4] R. Conte, Universal invariance properties of Painleve analysis and Backlund transfor- 
mation in nonlinear partial differential equations, Phys. Lett. A 134 (1988) 100-104. 

[5] R. Conte, Invariant Painleve analysis of partial differential equations, Phys. Lett. A 
140 (1989) 383-390. 

[6] R. Conte and K.-w. Chow, Doubly periodic waves of a discrete nonlinear Schrodinger 
system with saturable nonlinearity, J. Nonlinear Math. Phys. ? (2008) accepted. 
|http: / / arxiv.org/abs/0812 . 1 196| 

[7] R. Conte, A. P. Fordy, and A. Pickering, A perturbative Painleve approach to non- 
linear differential equations, Physica D 69 (1993) 33-58. 

[8] R. Conte and M. Musette, Painleve analysis and Backlund transformation in the 
Kuramoto-Sivashinsky equation, J. Phys. A 22 (1989) 169-177. 

[9] R. Conte and M. Musette, Link between solitary waves and projective Riccati equa- 
tions, J. Phys. A 25 (1992) 5609-5623. 

[10] R. Conte and M. Musette, Analytic expressio ns of hy drothermal waves, Reports on 
mathematical physics 46 (2000) 77-88. |nlin.SI/0009~022| 

[11] R. Conte and M. Musette, The Painleve handbook, xxiv+256 pages (Springer, Berlin, 
2008). |http://www.springer.com/physics/book/978- 1-4020-8490-4| 

[12] R. Conte and Ng T-w, Meromorphic solutions of a third order nonlinear differential 
equation, to be submitted (2009). 

[13] G. Darboux, Sur les equations aux derivees partielles, C. R. Acad. Sc. Paris 96 (1883) 
766-769. 

[14] J.D. Fournier, E.A. Spiegel and O. Thual, Meromorphic integrals of two nonintegrable 
systems, Nonlinear dynamics, 366-373, ed. G. Turchetti (World Scientific, Singapore, 
1989). 

[15] Mark van Hoeij, package "algcurves", Maple V (1997). 
http://www.math.fsu.edu/~hoeij/algcurves.html 



[16] A.N.W. Hone, Non-existence of elliptic travelling wave solutions of the complex 
Ginzburg-Landau equation, Physica D 205 (2005) 292-306. 

[17] A. Jeffrey and Xu S., Travelling wave solutions to certain non-linear evolution equa- 
tions, Int. J. Non-Linear Mechanics 24 (1989) 425-429. 

[18] A. Khare, K.0. Rasmussen, M.R. Samuelsen and A. Saxena, Exact solutions of the 
saturable discrete nonlinear Schrodinger equation, J. Phys. A 38 (2005) 807-814. 

[19] N. A. Kudryashov, Exact soliton solutions of the generalized evolution equation of 
wave dynamics, Prikladnaia Matematika i Mekhanika 52 (1988) 465-470 [English : 
Journal of applied mathematics and mechanics 52 (1988) 361-365] 

[20] N. A. Kudryashov, Exact solutions of a generalized equation of Ginzburg-Landau, 
Matematicheskoye modelirovanie 1 (1989) 151-158. [English : none?]. 

[21] Y. Kuramoto and T. Tsuzuki, Persistent propagation of concentration waves in dis- 
sipative media far from thermal equilibrium, Prog. Theor. Phys. 55 (1976) 356-369. 

[22] P. Manneville, Dissipative structures and weak turbulence (Academic Press, Boston, 
1990). French adaptation: Structures dissipatives, chaos et turbulence (Alea-Saclay, 
Gif-sur-Yvette, 1991). 



16 



[23] M. Musette and R. Conte, Analytic solitary waves of nonintegrable equations, Physica 
D 181 (2003) 70-79. |http: / / arXiv.org/abs /nlin.PS /030205l| 

[24] A. Pickering, A new truncation in Painleve analysis, J. Phys. A 26 (1993) 4395-4405. 

[25] W. van Saarloos, Front propagation into unstable states, Physics reports 386 29-222 
(2003). 

[26] H. Segur, Solitons and the inverse scattering transform, Topics in ocean physics, 
235-277, eds. A. R. Osborne and P. Malanotte Rizzoli (North-Holland publishing 
co., Amsterdam, 1982). 

[27] S. Yu. Vernov, Elliptic solutions of the quintic complex one- dimensional Ginzburg- 
Landau eq uation, J. Phys. A 40 9833-9844 (2007). 
|http: / /arXiv.org/abs /nlin.PS /060'20"60| 

[28] J. Weiss, M. Tabor, and G. Carnevale, The Painleve property for partial differential 
equations, J. Math. Phys. 24 (1983) 522-526. 



17 



